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Abstract 

°° (— )n 1 

Infinite series of the type ^ - 2 -^— 2-^1 ( — n, b; 7; y) are investigated. Closed-form 

n=i n ' 

sums are obtained for a a positive integer a = 1, 2, 3, ... . The limiting case of b —* 00, 

OO ,a_ \ 

after y is replaced with x 2 /b, leads to ^ -^-^| n, 7, x 2 ). This type of series 

71=1 

appears in the first-order perturbation correction for the wavefunction of the general- 
ized spiked harmonic oscillator Hamiltonian H = — 4-? + Bx 2 + -4r + "4r < x < 
00, a, A > 0, A > 0. These results have immediate applications to perturbation se- 
ries for the energy and wave function of the spiked harmonic oscillator Hamiltonian 
H = --£ 7 + Bx 2 + ± 0<x<oo, a,X>0. 

PACS 03.65.Ge 
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1 Introduction 

Recently, Hall et al |L], 0, [| have investigated some infinite series of the type 



oo 



E^Vit-W 2 ), 7 >"a = 2,4,... (1) 
' n n! 2 

n=l 

where i-Fi stands for the (Kummer's) confluent hypergeometric function defined in terms of 
associated Laguerre polynomials by 

n ( \ k I 

^ (7 + l) fe k\ (7 + l) n 
and (j) n , the shifted factorial (or Pochhammer symbol), defined by 

( 7 )o = l, ( 7 ) w = 7 ( T + l)( T + 2)...( T + n-l)= y J , n = l,2,.... (3) 

where 

if < < n 

-/')*■ = < (B -* ) " ... " " ( 4 ) 
0, if fc > n. 

Beside its mathematical relevance as generating functions for Laguerre polynomials [§], [|, 
the infinite series ([I]) appears in attempts to develop wavefunction perturbation expansions 
for the generalized spiked harmonic oscillator Hamiltonian || [7], || 

A d 2 9 A A 

H = H Q -\ = — r^+x 2 + — -\ < x < 00, A>0,a,A>0, (5) 

x a dx l x l x a 

which, in turn, has immediate applications in the study of the well-known spiked harmonic 
oscillator 0, |, § 

d 2 A 

H = — — - + x 2 H 0<x<oo, a, A > 0. (6) 

aar x a 

In particular, it helps to explain the abnormal behavior of the standard, weak coupling, 
perturbation theory || for Hamiltonians @. In the usual elementary sense the potential 
l/x a is singular whenever a > 0. However, the value a — | marks the division of the problem 
into two quite distinct regimes: we shall adopt the convention (following Harrell |TU| ) that 
the term "singular" means a > |; some writers such as Klauder |TT| use instead the term 
"supersingular" for this regime. Hall et al found || 0, ^| that for the Hamiltonian (H) 
defined on the one-dimensional space (0 < x < 00) with eigenfunctions satisfying Dirichlet 
boundary conditions, that is to say, with wavefunctions vanishing at the boundaries, the 
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matrix elements of the operator x Q , with respect to the exact solutions of the Gol'dman 
and Krivchenkov Hamiltonian Hq, namely 



with exact eigenenergies 



1 



E n = 2(2n + 1 ), n = 0,1,2,..., 7 = 1 + -VT+IA 



are given explicitly by the following expressions 



x ™ ~ Wn~fW~ V ^!^r 3F2(_m ' 7 " i ' 1 " 2 ' 7 ' 1 ~ n ~ 2 ' 1} (9) 

and they are valid for all values of the parameters 7 and a such that 7 > |. Furthermore, 
the matrix elements of the Hamiltonian (^) are given by 

H mn = 2{2n + 7)^+(-l) m+ "AA/ (7) " (7 , )m F( J " $ ( f )n x 

V n\m\ (7) n r(7) 

„ . a a ol , 

3 F 2 (-m, 7 --,l--;7,l- 2 - ^1)- (10) 

where 6 mn denotes the Kronecker delta. Of particular interest are the elements 



It is known that the first correction to the wavefunction by means of standard perturbation 
techniques leads to 



n=l 



where ip n (x) and H 0n are given by (|?p and (JTT|) , respectively. Thus, the first correction to 
the wavefunction of the Hamiltonian fl5|) is given by 

,(1)/ \ 1 r (7-f) 7-1- 4\^(f) nl r? I 2\ / 10 \ 

ipn ix) — -= . x' 2 e 2 \ -Jl - X \F\{— n, 7, x ). (12) 

2 v / 2r( 7 ) v T T (77 ^ r, ^ L u ; 



n=l 



n n\ 



Very little indeed is known about the infinite series which appears in (12) for arbitrary 



values of a. The first attempt to sum this series, for the very particular case a = 2, is due to 
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Toscano who proved, by means of extensive use of calculus of finite difference [13[], and 
in terms of associated Laguerre polynomials (|J), that 

(n-l) 



E 

n=l 



(7)» 



1 4 7 " 1) (l/)=^(7)-logl/, 7>1 



(13) 



where VKt) is the Digamma function defined by VKt) = dlog j r ^ 7 ^ . Another attempt to sum 
this series was due to Aguilera-Navarro and Guardiola ||, who reported that they met some 
serious convergence difficulties even in the case a = 2 when trying to get the first correction 
to the wave function. Recently, however, Hall et al gave a systemic method for summing this 
series in the case of even a = 2, 4, 6, ... . They proved for ^ = 2 + m, m = 0, 1, 2, . . . that 

(2 + m), 



E 



n=l 



n n\ 



L Fi(-n; 7; x ) = ^(7) - logs 



(m 



+ DE 



fc=0 



(fc + 1) 2 



( 7 - 1) r7 -2- 



X ' 



where -^^(y) defined by (fj), and furthermore, in the case of 



00 ^ 

V" -iFi(-n; 7; x 2 ) = ^(7) - logx 2 



(14) 



(15) 



n=l 



In this article, we continue our search for closed-form sums of the perturbation series appear- 
ing in (12) by developing more effective methods for the even case and by introducing new 
methods for the odd integers a = 1, 3, 5, ... . We achieve this by studying a more general 
class of infinite series involving the classical hypergeometric functions 2F1 of Gauss, namely 



z — ' n n< 



n=l 



a 

7 > — , a 
1 2' 



1,2,... 



This leads, as a result of the limiting case of, 



x 



lim 2 F 1 (-n,b;r, , 

0— >oo 



L Fi(-n;7;x 



(16) 



(17) 



to our series in (|T|). The family of infinite series (|16|) and (P are related to the theory of 
generating functions as well to the theory of the inverse Laplace transform: as a side goal we 
shall explore these links too. The paper is organized as follows: in Sec. II we demonstrate 
that the infinite series (|16D converges for all y > and 7 > |. Furthermore, we prove that 
the series ([16]) has the following integral representation valid for 7 > b 



00 (a\ 

E^?^i(-n,6;7;y) 







r(7) 



n=l 



n n\ 



2 r(6)r(7 - b) 



x 



ft 



t*-i(l _ _ ^) 3 F 2 (1 + - 1, 1; 2, 2; 1 - yt) dt 
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In order to evaluate the integral just mentioned, we develop, in Sec. Ill, analytic expressions 
of the series 3-F 2 (l + ~, 1, 1; 2, 2; z) for arbitrary values of a and \z\ < 1 and we shall prove 
that 



z 3 F 2 (a+l,l,l;2,2;z) = ^J £ J. dt, a ^ 0. 

Sec. IV is devoted to the even case a = 2, 4, 6, . . . , m + 2, m = 0,1,2,... where we shall 
prove, for y > 0, the general formula 

°° (m + 2) 

V "^(-n, 6; 7; 2/) = ^(7) - 4>(b) - log(y) 



n=l 



,_„ (2)*(2)* 



1 ^(-A;, 2 - 7; 2 - 6; -) - 2 F 1 (-A; ) 1 - 7; 1 - 6; 1 



(6-1)2/ ?/ ?/ 

which leads, as we shall demonstrate, to ( |i~4"D as a special limiting case. In Sec. V, we shall 
find closed-form sums for the series ( Jl6| ) for the cases a = 1,3,5,..., and hence we shall 
prove that 

^2*i(-n, b; 7; y) = 2 log(2) + -y 3 F 2 (1, 1, 1 + 6; 2, 1 + 7; y) 

n=l 

2r( 7 )r(6 + i) 1131 
~ r(6)r( 7 + i)^ sF2(1 ' 2' 2 + 6; 2' 2 + 7; y) 

which generate closed form sums for arbitrary odd-integer values a = 1, 3, 5, ... . Closed-form 
sums for the series (||) then follow as a special limiting cases of the result just mentioned. 
Finally, in Sec. VI we study some applications to the wave function perturbation expansion of 
the generalized spiked harmonic Hamiltonian: we shall give a complete first-order correction 
to the wave function in the case of the radial harmonic oscillator as well as for the spiked 
harmonic oscillator Hamiltonians. Some other applications to the theory of inverse Laplace 
transform and generating functions will also be discussed. 

It is clear that the functions 3 F 2 , and 2 i*i , mentioned above, are special cases of the 
generalized hypergeometric function 

p 

pF q (ai,a 2 ,..., Op] fit, fa,..., p q ;z) = jj, 



k\ 

Pj)k 



This series converges for all z if p < q + 1, and for \z\ < 1 i£ p — q + 1. It has been assumed 
that all parameters have real values except for the /3j,j = 1,2, ... ,q, none of which is equal 
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to zero or to a negative integer. If p = q + 1, the series is absolutely convergent on the unit 
circle \z\ = 1 if 



£>-Z> >0 - 



(19) 



i=l 



If one or more a% are negative integers, the series terminates as a result of (Q) and convergence 
does not enter the discussion. If p = q, the series (181) converges for finite z. 



2 Integral Representation and the Convergence Prob- 
lem 



In this section we prove the convergence of the infinite series (16) and develop an integral 
representation suitable for our calculations. It is important, however, to notice for y = 0, 
the infinite series (0) for a > 2 indeed diverges. Since 2-^i( — n , b; 7; 0) = 1 and (f) n +i = 
f (f + 1)„, we have 



(20) 



= = > 2 d + ?, 1, 1; 2, 2; i) = m - * (i a 



nw 2^ (2) n (2)„n! 

n=l n=0 



2 U v 2' ' ' ' ' ' TK ' Ty 2 
which is absolutely convergent for a < 2. Furthermore, for y = 1, it is known that 

(7 - b) n 



2 F 1 (-n,6;7; 1) 



(7), 



as a result of Chu-Vandermonde theorem ||15|| . Thus 



(f )n (7 - b) n 1 = q( 7 - 6) A (1 + f )n(7 - & + 1 



^-f n (7) n n! 

n=l v ' 



7 



n=0 



(7 + l) n (2) n (2) n n\ 



^— ^ 4 F 3 (1 + " 7 - 6 + 1, l, l; 7 + l, 2 , 2; 1) 
7 2 



(21) 



which is absolutely convergent for | < 1 + b due to fll9f). For arbitrary y > 0, and in order to 
sum the infinite series (fT6j) for 27 > a > 0, we require a suitable integral representation of the 
hypergeometric function 2^i( — n, b; 7; y) over an appropriate contour, in order to interchange 
summation with integration and thereby readily conclude the absolute convergence of the 



series just mentioned. We find the Euler's integral representation |16 

1 



2*1(0, 6; 7; y) 



r(7) 



r(6)r( 7 - b) 



t b -\i-ty- b - 1 (i-yt)- a dt 



\7-6-l 
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under the conditions Re(j) > Re(b) > 0, and \y\ < 1 to be most advantageous for achieving 
this end. For a = —n, these conditions can be relaxed as a result of (Hf), we have 



i 



2 F 1 (-n,b;r,y) = T ^^Jt b -\l-tr- b -\l-ytrM (22) 

o 

for arbitrary value of y > 0, which, when substituted into the summation of (|16l), yields 

l 

77. — 1 ft — 1 q 

The evaluation of the this last infinite sum, involving integrations over the open interval 
(0, 1), is achieved by examining the summation of the integrand, namely 

^V^i - ty-»-\i - ytr = t»-\i - ty- b -> y: - yt) n 

I'll. till. 

and demonstrating that it has an Li(0, l)-majorant. Hence, the existence of such a majo- 
rant shall permit us to interchange summation with integration, as result of the Lebesgue 
Dominated Convergence Theorem. To arrive at such a majorant, we continue by noting that 

Oi OL 

= -(1 - yt) 3 F 2 (l + -, 1, 1; 2, 2; 1 - yi) (24) 
which is convergent for 1 1 — | < 1 as a result of (HI). Thus 

I £ il22f^i(i _ t )T-fc-i(i _ yt y\ = ^t b -\l - ty- b -'\l - yt\ | 3 F 2 (1 + 1, 1; 2, 2; 1 - yt)| 
z — ' n. n! z 2 

< ^(a^-^l-t) 7 - 6 - 1 , te(0,l) (25) 



n=l 



where the absolute convergence of 3F2 and also |1 — yt\ < 1 were made use of. The most 
important aspect of inequality ( p5|) is the appearance of the Li(0, l)-function £ 6_1 (1 — t) 7-6-1 , 
for 7 > b, of variable t majorizing the series 

Y k4| t 6-l(l_ t )7-6-l(l_ y£ )nl 

z — ' n n! 

n=l 
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This aspect justifies the evaluation of summation (^) by means of the Lebesgue Dominated 
Convergence Theorem. Thus, we have specifically, for 7 > b, 

E(f )« n / 7 \ « T(7) 
>-| ^^»)'i rt W7 -n '' 

1 

£ 6-i (1 _ t) 7-6-i(! _ yt ) 3 F 2 (l + |, 1, 1; 2, 2; 1 - yt) dt (26) 

For the computation of this expression, we need analytic expressions for the hypergeometric 
functions 3-F 2 (l + |, 1, 1; 2, 2; 1 — yt) : these are obtained in the next section. 



3 Analytic expressions for 3^2(1 + § , 1, 1; 2, 2; z) 



In this section we present some introductory results concerning the analytic sum of 3^(1 + 
1, 1, 1; 2, 2; z) for later use. We start first with an interesting relation between a Gauss 
function and an elementary function, not usually shown in specialized tables. 

Lemma 1. For \z\ < 1 

jpf i 1 i o ^ j-1og(l-z), ifa = 
z 2 F 1 [a + 1,1; 2; z) = < , : (27) 

Proof. From the series representation of the the hypergeometric function 2 F\, (0), we have 
2 F 1 (a + l,l;2;z) 



However, from the identity 

dz 



oo 

=£ 

n=0 


(a + l) n (l)„s n 


(2)n 


n! 


oo 

= E 

n=0 


(fl+l)„^ 

n + 1 n! 




oo 

= E 

n=0 


z n 

(a + l) n — r 
n! 


Y> n(a + l) n z n 
' n+1 n! 

71=1 


= 1^01 


> + l;-;z) 


2 ^ (3) B n! 


= i-Fol 


>+ l;-;z) 


(a + l)z 
- 2 2 F 1 (a + 2,2;3;z) 


a, b; c; 


ab „ 

c 


[a + 1, b + 1; c + 1; z) 
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we have 



2 F 1 {a + 1, 1; 2; z) = x Fo(a + 1; -; 2) - z^- 2 F x {a + 1, 1; 2; *) 

u2 



which yields 



d 

dz 



z 2 FAa + 1,1:2; z) 



(1-4 



-(0+1) 



where we have used the identity i-F (a + 1; — ; z) = (1 — z) ( a+1 ), |z| < 1. Therefore, 

2 Fi(a + l,l;2;z) = / (1 - t)~^dt 
Jo 



z 2 j 



For a = 0, we have z 2-^1(1, 1; 2; z) = — log(l — z), while for a 7^ the result 

z 2 F l (a + l,l;2;z) = — + ' 



a a(l — z) a 



follows immediately by elementary integration. 



□ 



Lemma 2. For \z\ < 1 

z 3 F 2 (a4 1-1.1:2.2:: 



Li 2 (z) 



if a = 



where Li 2 (z) is the dilogarithmic function defined by 



(28) 



Proof. From the series representation fll8|) of the hypergeometric function 3F2, we have 

(o+l) n (l) n (l) n 2» 



,F 2 (a + 1,1,1; 2,2;*) = ^ 

n=0 

oo 

=E 

n=0 

00 

= E 



(2) n (2)„ n! 



(a + l) n (l) ra z" 
(n + l)(2) n n! 

(a + l)„(l) n z" 



-E 



n(o + !)„(!)„ 2 P 



n=0 



(2)n n\ ^ (n+l)(2) n n! 



2 Fx(a + 1, 1; 2; z) - EL±^ 3 F 2 (a + 2, 2, 2; 3, 3; z) 
2 Fx(a + 1, 1; 2; z) - ^3*2(0 + 1, 1, 1| 2, 2; z) 
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which leads to 



d 
dz 



z 3 F 2 (a + l,l,l;2,2;z 



^(a + M^z,, 



or 



z 3 F 2 (a + 1,1,1; 2,2; z) = [ 2 F l (a + 1, 1; 2; t) dt. 

Jo 

For a = 0, and using Lemma 1, we get 

r* r* loe-d - 1) 

z 3 F 2 (l,l,l;2,2;z) = / 2 F 1 (1, 1; 2; t) dt = - / ^dt = U 2 (z) 

Jo Jo J 

while for a 7^ 0, we have by means of Lemma 1 that 



z 3 F 2 (a + 1, 1, 1; 2, 2; * 



1 y 1 -* r - 1 



a t a (l -t) 



dt 



as required 



□ 



The definite integral in Lemma 2 is easy to evaluate. It follows immediately, for example, in 
case of a = 1 that 



z 3 F 2 (2 } 1, 1; 2, 2; z) = z 2 F X (1, l;2;z) = -J jdt = - log(l - z) (29) 
as expected. For a = ~, 

z 3 F 2 (|, 1, 1; 2, 2; z) = 4 log 2 - 4 log(l + v 7 !^), (30) 

for a = 2, 

1 r 1 

(31) 



z 3 F 2 (3, 1, 1; 2, 2; z) = - \ log(l - z). 



2 1 - z 2 

Similar analytic expressions can be obtain for arbitrary values of a. However, for a = 
m + 2, m = 0,l,2,..., the following lemma which can be regarded as a limiting case (as 
b — > 1) of Luke's reduction formula |U 



z 3 F 2 (a,b, l;c,2;z) 



(c-1) 



(o-l)(6-l) 



2 Fi(a-l,6-l;c-l;z) - 1 



Izl < 1, 



can be used. 
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Lemma 3. For c ^ 1, and I^M < 1, 



z 3-^2(0 + 1,1, l;c, 2;z) 



(c-1) 



fc — a — 1"! 



c-1 Vz-l 



\ Z 

3 F 2 (c-o,l,l;c,2; 



z- 1' 



Proof. For the proof of this lemma see Lemma 3 in [|T] 



log(l - z) 
(32) 

□ 



Although, the right hand side of ( j32T) still in terms of 3F2, it is indeed more straightforward 
to deal with in the case of a — m + 2, m — 0, 1, 2, ... Indeed for c = 2, we have 



2 3 F 2 (m + 3,l,l;2,2;z 



m + 2 



(m + r 



3 F 2 (-m,l,l;2,2; 



z — 1 J z — 1 

where 3_F 2 (— m, 1, 1; 2, 2; is now a polynomial of degree m, therefore 

1 



bed - z) 



z 3 F 2 (m + 3,l,l;2,2;z) 



m + 2 



m + 1 



z — 



1 £ 



ra) n (l) n (l) n 



n=0 



(2) n (2)„n! V^-l 



- logfl - z) 



(33) 



which leads, for example, to 



z 3 F 2 (4,l,l;2,2;z) 



1 / z 



2z 



2\l-z 1 - z 



- log(l - 2) 



The case of unit argument z — 1, the hypergeometric function 3 F 2 (a + 1, 1, 1; 2, 2; 1) diverges 
for a > 1, while it converges absolutely for a < 1. With the strict inequality |z| < 1, the 
function converges absolutely for arbitrary values of a. 

4 Closed-form sums for the even case a = 2, 4, 6, . . . 

In this section we develope closed form sums for the series ([16]) for a = 2,4, 6, . . . which 
generalize the methods and results of Hall et al @, [|. The odd case a = 1, 3, 5, . . . will 
be discussed in the next section. 

Lemma 4. For 7 > b and y > 

00 j 

V - 2 Fx(-n, 6; 7; 2/) = - " log?/ (34) 
z — ' n 



n=l 
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Proof. From a = 2, using fl2"9|), the integral representation fl26|) leads to 



53 -2^1 (-n, 6; 7; 2/) 



r( 7 ) 



n=l 



r(6)r( 7 - 6) 



1 

y t 6 -^! - 1)^- 6 - 1 iog(t)dt 



r(7) 



r(6)r( 7 



1 

-^iog(y) J t h -\i-ty- h - l dt 



The second integral on the right hand side follows immediately from ( p2] ) by setting n = 
and leads to 



1 



r(fe)r( 7 - b) 
r( 7 ) 



(35) 



while the first integral can be computed by differentiating (|35|) with respect to 6, which leads 

1 1 
to the difference of / ^(l - f)7-*>-i log(t) dt and / ^(l - i)"?- 6 " 1 log(l - t) dt, where the 


latter can be regarded as the derivative of (|35| ) with respect to 7. This yields 



dfe V r (7) 



c?7 V T(7) 



which leads, after some simplification, to fl34| ) as required 



□ 



It is clear from fl34| ) that the condition 7 > b can be removed by analytic continuation 
consequentely the special case of the infinite series fll^D follows by replacing y with ^- in 
and then taking the limit as b approaches 00, this yields 



V"-iFi(-n;7;x 2 ) = lim 



n=l 



^(7) — ^(b) + log 6 — log a; 2 



as a result of ([[]])■ However, from the integral represention of the Digamma function 

00 

^(6) = k)g(6) + y [r 1 - (1 - e^) _1 ]e" W dt, »(&) > 



it follows that lim (^(6) — log 6) = and (O) follows immediately. Similar results for lemma 

(4) can be carried out for the values of a = 4, 6, . . . using the integral representation ( ^6l) and 
the analytic expressions for 3F2 as developed in Sec. 3. However, we shall adopt a different 
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approach which allows us to obtain a general formula for the sum of the series (|T6|). First, 
we generalize Buchholz's identity for Laguerre polynomials |T 



2^ r^^^^l) L " ^ = 2_, ZTn77~T~T\ 7 + 1; y) = y , (36) 



r(n + 7 , ... , 

n=0 V ' ' n=0 



n!r( 7 + l) 

for 7 + z/> -1, i/^ 0,1,2,.... 

Lemma 5. For 7 + z/>0, o + z/>0, and y > 



E 

n=0 



— 2 iq(-n,6; 7 ;y) = W)N — , _^ 



r(6) r(!/ + 7)' 



(37) 



Proof. From the Laplace transform representation [|18|] of 1 iFi(— n; 7; yt) with unit argu- 
ment, we have 

00 

2 Fi(-ti,6;7;j/) = fj^ J e'^ 1 xFi{-n; 7; yt) dt, 



Thus 



E 

ra=0 



n! 



- 2 Fi(-n,&;7;y) 



r(6) 



r(7) 



£ 

n=0 

00 



-iFx(-n;7; yt) 



fit 



r(6)r( 7 + z/ 

r( 7 )r(6 + */ 
r(6)r( 7 + z/ 



■y" y e"*^- 1 dt, 7 + 1/ > 


•y", b + u>0 



(3J 



where (J38T) follows from (|36|) and the last equality follows from the definition of the Gamma 
function. □ 

With the availability of Lemma 4 and Lemma 5, it is relatively easy, using (E^), to find the 
closed form sums for the series (|T6|) and consequentely the infinite series (|l|) for a = 4,6, ... . 
Indeed, for a = 4, we note that (%) n = (2)„ = (1 + n)(l) n , thus 

V ^y 2 i ? i(-n, 6; 7; 2/) = £ ^(-rc, 6; 7 ; y) + V 2 i ? i(-n ) 6; 7; y) 



n=l 



n=l 



n=l 



The first series on the right hand side is summable by means of Lemma 4, while for the 
second series, it is enough to take v — — 1 in the Lemma 5. This leads to 



00 (2) — 1 

£ r3T2^i(-^, 6; 7; 2/) = ^(7) - - lo S2/ + 



71=1 



n n; 



(6-i)y 



1, 6>1,7>1 
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and consequently 



jt ^-M-n; 7 ; * 2 ) = ^(7) - log* 2 + ^ - 1 



71=1 



follows by means of fll7| ) and lim6_ i . oo ('0(6) — log(6)) = as mentioned earlier. Similarly, for 

a = Q and (3) n = |(2 + 3ri + n 2 )(l) n , we have 

OO ,_N OO _ OO 

£ ^a^iC-n, 6; r,v) =Y, - 2 *i(-n, 6; 7; y) + £ £ a fi(-n, 6; 7; v) 

71=1 
^ OO 

+ 9 E n 2-Pi(-n,6;7;y) 



n=l 



n=l 



71=1 



Closed form sums for the first and second infinite series on the right-hand side can be found 
as mentioned above, for the third infinite series it is enough to take v = —2 in Lemma 5. 
These results lead to 

£ %^(-".»; r, v) = m - m - io g!/ + + j^l - 3 - (39) 



n=l 



2(6-l)(6-2)j/2 ( 6 _i) y 2 



for 7, 6 > 2, and consequentely 

OO 

£ ^T 2 ^ 7; = ^(7) - log x 2 + 



(3)n 2 , , , v , 2 , 1 (7-l)(7- 2) ,7-1 3 



+ 



71=1 



n n! 



ar 



(40) 



by means of (|1~?D . For arbitrary values of | = m + 2,m = 0,l,2,..., a general expressions 
for the closed form sum of the infinite series (|16|) can be obtained by means of the following 
Lemma. 

Lemma 6. For 7 > m + 2 and b > m + 2, where m = 0, 1, 2, . . . 



n=l 



fm + 2) 



n n\ 



- 2 Fi(-7i, 6; 7; y) = ^(7) - - lo s(y) 



+ ( m +l)]T 



fc=0 



-m) fc (l)fc 
(2) fc (2) fc 



;7 }) *F x (-k, 2 - 7; 2 - 6; -) - 2 F 1 (-A;, 1 - 7; 1 - 6; i) 

(o-i)j/ y y 



(41) 



Proof. From the integral representation (|26[) and the analytic expression (^), we have 



E 

71=1 



[m + 2), 



n n! 



2 Fi(-n,6;7;y) 



r(7) 



r(6)r( 7 - 6) 



t 6 - 1 (l-t)^- 6 - 1 log(yt) rft 



(m + l)r( 7 ) 

r(6)r( 7 - 6) 



(2) k (2) k k\ J y ' \yt )\ ytj 
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The first integral on the right hand side is easy to evaluate by a similar proof to Lemma 1 
and yields 

~ r(6)r(7- b) jf " ty ^ logiyt) dt = Hl) ~ m ~ logiy) 

The second integral, which we denote it by I m , can be written as (since (1)„, = n\) 



r(i)r(7-i)g (2)*(2)» yj ' V V* 

(m + l)r( 7 )^(-m) fc (l) fc /* i 



r(6)r( 7 -6)^ (2) fe (2) fc 



Since 



1 yt> ^ l\ K yt> 
y i=o y 



» i — /,• i , i 

finite number of terms, 



we have 



_ (m + l)r( 7 ) ^ (-m) fc (l) fc l A (-k)i 1 Z" 1 ft-j-a/, _ -w-6-1 

im " r(6)r( 7 - 6) ^ (2) fc (2) fc /! Wo 1 j 

(™ + l)r( T ) f (-m) fc (l) fc A (-*), 1 f 1 



\- (-m)fc(l) fc \- (-Apt 1 /" , 6 -i-l/ 1 _,w-6-l 
£± (2) fc (2) fe 2. «l yl J Q t 1 ^ 



r(6)r( 7 -6)^ (2) fc (2) fe ^ /! y« 

where each of the integral on the right hand side can be evaluate by means of Beta function, 
this leads to 

(m + l)r( 7 ) (-m) fe (l) fc 1 A (-fc); 1 r(6 - / - i)r( 7 - b) 

m r(6)r( 7 -6)^ (2) fc (2) fe y^ u y l r( 7 -J-i) 

(m + l)r( 7 ) (-m) fc (l) fe A (-*;), 1 T(6 - fc)r( 7 - 6) 



r(6)r( 7 -6)^ (2)*(2)* ^ Z! y fc r( 7 -*) 
Using the identities 

(-i) ; r(a) (-iyr(a-i) 

1 (a — 1) = — : — , l(a — I — 1) = z : 

1 ' (l-a)i (2 -a), 

we obtain by means of the series representaion (|20|) of 2 -Pi that 

(m+l)( 7 -l) y- 2 - 7 ;2-6;i) 

. (— m)*.(l)i. . 1 , 

-^E^^t-M-r.i-^) 

which complete the proof of the Lemma. □ 
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The limit cases 



and 



lim 2 F x {-k, 2 - 7; 2 - 6; — ) = 2 F (-k, 2 - 7; -; -— ) 



6— >oo 



6— >oo 

which follows as results of 

T(2 - 6)6 



lim 2 F 1 (-k, 1 - 7; 1 - 6; — ) = 2 F (-£;, 1 - 7; 



ar 



lim 



b^ooT(2-b + n) v 
respectively, lead us to conclude that 
" (m + 2) n 



-l) n , and lim 



T(l - b)b n 



b-*oo r(l - 6 + n) 



-iy 



n=l 



n n! 



lFi(-ti; 7; x ) = ^(7) - log(a: ) 



+ (m+l)^ 



fc=0 



^-m)fc(l) fe 
(2) fc (2) fc 



' ' " 1} 2F (-k, 2 - 7; -3) - 2F (-k, 1 - 7; 



ar 



.2;" 



which agrees with ([16]) through the identity |I5 



4 a " n) (2/) = LJ i L 2F (-n, -a; -; — ), for n = 0, 1, 2, 
n! y 



(42) 



(43) 



5 Closed- form sums for the odd case a = 1, 3, 5, . . . 



We now turn our attention to the case of a = 1, 3, 5, ... , namely | = m + |, m = 0, 1, 2, 

In this case, we don't have, indeed, a concise formula similar to (0), and our results will 
depend on a recursive approach. For example, in the case of a = 3, we notice that (|) n = 
(§)n — (1 + 2n)(|) n and therefore 

> -^72^i(-n, 6; 7; y) = > j-^— 2 i ? i(-n, 6; 7; y) 



n n! 

n=l n=l 



n ni 



00 / 1 - 



V ^Fii-n, 6; 7; y) + 2 V 2 F x (-n, 6; 7; y) (44) 



n=l 



n=l 



nl 



The second series on the right-hand side is summable by means of Lemma 5 which leads, for 
the case of v — —5, to 



00 ^1 
nl 



n=l 



(45) 
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Similarly, in the case of a = 5, and since (-)„ = (1 + |n)(|) n , we have 

E %^i(-n, 6; 7;y) = E (1 + ^ )( i ) " 2 F 1 (-n ) 6; 7 ; y) 
z — ' n n! z — ' n n! 

n=l n=l 



E ^FxC-n, b ]T ,y) + iYl ^r^iC-n, 6; 7; y) (46) 



n n\ 3 ^— ' n! 

n=l n=l 



The second series on the right-hand side again summable by means of Lemma 5 which leads, 



in the case v = — |, to 



®» w( h \ r (fe-|) r( 7 ) i 1 (A7 , 

g— 8 F 1 (-n,6; 7 ;y)= r(fc) ^ _ |} — - 1, 7 >- (47) 



while the first series is summable by means of (fig), if of course, we have a closed form sum 
for the infinite series 



OO / i ' 



E^fU-^t;* 2 ) (48) 



n n 

n=l 



Equivalently, for the cases a = 7, 9, . . . which depend on closed-form sums for the infinite 
series (01). 



Lemma 7. For 7 > 6 



V ^■2F 1 (-n, b; 7 ; y) = 2 log(2) + -y 3 F 2 (1, 1, 1 + 6; 2, 1 + 7 ; y) 
' n n! 7 



n=l 



Proof. From the integral representation (|26|) and the analytic expression (^) we have using 

m 

E %2^(-n, 6; 7 ; y) = 2 log(2) - ^ ^(l " *) 7 ^ Ml + v^) dt 

^ n n\ T(6)r( 7 - 6) J 

The integral on the right hand side can be performed analytically and leads to 

J! tb " (1 - V' + VSD * = 4 r(i+K +7 ) ^ r(6 + + 7) 
3 F 2 (i, 1, 1 + 6; 2, 1 + r , y) - 2r(i + 6)r(i + 7 ) 3 F 2 (i, i i + 6; |, i + 7; a ) 

which the reader can verify using symbolic software such as Maple V. After some careful 
simplifications, this leads to the result of the lemma. □ 
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The limit cases 



x 



lim 3 F 2 (1, 1, 1 + 6; 2, 1 + 7; — ) = 2 F 2 (1, 1; 2, 1 + 7; ^ 

6->oo 



and 



i im r ^ +6) 



6^00 y/bT(b) 

leads us to conclude for a = 1 

00 / 1 - 



n/ 11 ,31 - 2 
3^2(1, -, - + 0; -, - + 7; — J = 2^2 



' 2' 2 ' 2' 2 ' " ~b ' ~ ^ ZVi ' 2' 2' 2 



13 1 2x 
F 2 (l, -;-,- + 7; x 2 ) 



V ^a^i(-n; 7; x 2 ) = 2 log(2) + - 2 F 2 (1, 1; 2, 1 + 7; x 



n=l 



7 

2xr( 7 ) 



^2(1,^ i + 7;x 2 ), for 7 > J. (50) 



r(7 + i)^"""'2'2'2 

This can be used now to obtain a closed form expressions for the infinite series ([!]) for 
a = 3, 5, ... , indeed we may now conclude using ( fi~5|) and fl47|) , that for a = 3 



V 7; x 2 ) = 2 log(2) + - 2 F 2 (1, 1; 2, 1 + 7; * 2 ) 



n=l 

Similarly, for a = 5 

'n n! 



7 



2aT(7) 



1 3 1 



2r( 7 ) 1 



r(7+F F2(1 '2 ; 2'2 + 7; ^ + r(7-i), 



V^ii^-n; 7; x 2 ) = 2 log(2) + - 2 F 2 (1, 1; 2, 1 + 7; x 



2, for 7 > 



(51) 



n=l 



2T( 7 ) 1 2T( 7 ) 1 



2 ^ r (7) pn 1.3 1 , 1 2T(7) 1 

T(7 + |) 1 ' 2'2'2 +7,X j + r( 7 -|)x + 3r( 7 -|) x* 3 

Corresponding expressions follows in a similar fashion for a = 7, 9, ... . 



, for 7 > -. (52) 



6 Applications 

In this section, we shall discuss some applications for the results developed in Sec. 4 and 
Sec. 5. For the Hamiltonian (|5]) and @, it is known for the case of a < |, the potentials 
are non-singular, and one should expect that the standard perturbation theory could be 
applied without any difficulties. This is the case for the eigenvalue expansions, where the 
resulting expansion for the energy is easily obtained and agrees with numerically determined 
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eigenvalues for small A. Indeed, It has been shown by Hall et al [0] that the eigenvalue 
expansions, for the ground state energy, are given by 



Eo(a,y) = 2 7 + 



— r2( l P \ 2 4 F 3 (1, 1, - + 1, - + 1; 7 + 1, 2, 2; 1) + . . . 

I67 T 2 (7) 4 ^ ; ) 2 t , 2 -r , / t , , , , t 



r(7) 



(53) 



and valid for all a < 7 + 1, where 7 = 1 + Wl + 4A In particular, for a = 1 and 7 > \ 



^(1,7) = 2 7 + ^y^A 



3 3 



^^f^(l,l,§,§ + l;7 + l,2,2;l) + 



which leads for the spiked hamonic oscillator Hamiltonian, i.e. 7 = | or A = 0, to 
E (l, 7) = 3 + 1.128 379 167A - 0.077 890 972A 2 + . . . 



The eigenfunction expansions, however, cannot be always obtained with the same ease, nor 
do they have a similar form. To illustrate this point, the first-order expansions of perturbed 
wavefunctions, for x > 0, can be obtained using (|12"D and the results of Sec. 4 and Sec. 5 by 
the following expressions: for a = 1, 



^\x) 



1 r( 7 -|) i _j 

- -x y 2 e 2 



2V2r( 7 ) v / T(7)' 



x 



21og(2) + - 2 F 2 (1,1; 2,1 + 7; x 2 ) 



7 



2xr( 7 ) 
r(7TI 



- 2 F 2 fl, 



13 1 

2' 2' 2 



+ 7 ;x 2 ) 



1 



for 7 > -. (54) 



For a = 2 



^ 1} (x) 



V / 2(7-l) V / T(7) 



a; 2 e 2 



log 2 - -ip(j) 



, for 7 > 1. 



For a = 3 



1 r( 7 - 



x 7 2 e 2 



2v / 2r(7) v / T(7) 
2^(7) Hi 
r( 7 +±) 2 21 '2'2'2 



x 



21og(2) + - 2 F 2 (1,1; 2,1 + 7; x 2 ) 

7 



2r( 7 ) 1 



^- 2 F 2 (l, -;-,- + 7; x 2 ) + 2 



for 7 > -. 

' 2 



For a = 4 

#(x) 



x ' 2 e 2 



2V2( 7 -2)( 7 -l)v/f(7) 



log(x 2 ) - ^(7) 



7-1 



+ 1 



.1- 



(55) 



(56) 



for 7 > 2. (57) 



and similarly for a = 5, 6, 7, ... . In these expressions, the coefficient of the exponential 
term is no longer a polynomial in x. These are quite surprising results and demonstrate 
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clearly that the forms of the perturbation series for eigenvalues and wave functions are very 
different. 

It is clear with A = 1(1 + 1), where I — 0, 1,2, ... is a quantum number, the Hamil- 
tonian (||) can be regarded as the perturbation of the radial harmonic oscillator, while for 
A = (l + ~(N — l))(l + ~(N — 3)) it can be regarded as perturbation of the harmonic oscillator 
in N- dimensions, in either case (|53|)-(|5"6D provide the first-order perturbation correction to 
the wavefunctions for the cases of a = 1,2,3,.... There is another important point that 
follows from the properties of Pochhammer symbol (a + l) n = ^-^(a) n and Lemma 5 in Sec. 
4: closed-form sums for the infinite series (|i~6|) , and consequentely (H) for arbitrary a > 
and 1,2,3,... depend only on closed-form sums for series in the case of < a < 1. Al- 
though we know in this case the potential is non-singular and the Hamiltonian has accurate 
eigenvalue expansions that agree with the numerical results, the eigenfunction expansions, 
still however, have non-polynomial expansions. 

A side outcome of this work is that our results, especially (|2^), lead us to obtain some 
new formulas for the Laplace inverse transform that are not mentioned in the known tables 
dealing with hypergeometric functions [F2(J]. Hall et al |TJ introduced the following Laplace 
inverse transform representation of the infinite series ([[]) 

c+ioo 

£%iiM-n,7,* 2 ) = ^ / eV^l--) 3^2(1, l,l + ^2,2;l-^)d S . (58) 
' n n\ 2 2m Is 2 s 

n=l J . 

C—lOO 

The integral representation fl26|) allow us to conclude that the Laplace inverse transform of 
s -7(i _ ad) 3 p 2 (l, 1, 1 + f ; 2, 2; 1 - ^) with unit argument is 

£ ~ 1 ( s "< 1 - t» ^ 1 + -2 ■ x 2; 1 - t») = & mk^i x 

i 

j t b -\l - ty- b -\l - ^) 3 F 2 (l + |, 1, 1; 2, 2; 1 - ^) dt (59) 
o 

which leads, for example, in the case of a = 2, to the known formula 

^(^Qogs-logx 2 )) = ^^(7) -logx 2 ] (60) 

There are more interesting formulas for Laplace inverse transforms which are not usually 
mentioned in Tables: for 7 > | and x > 



£- 1 ^- 7 log(^ + ^ 



1 V'(7) . x 13 1 2 

2f^ + ft7T|) 2F2(1 '2 ; 2'2 +7;x) 

" 2r(7 + i) 2jF2(1;1;2 ' 1 + 7;x2) - (61) 
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Indeed, from (|58|), we have in the case of a = ~, using fl30D, that 

oo i c+ioo 

^^Ut,, 2 ),^ / e^(l--)aF 2 (?l,l;2,2;l--)d a 
n n\ 2 2m J s 2 s 



n=l 

c—zoo 



c+ioo 

^ ^ e s s^ ( 2 log 2 - 2 log(l + -^=) ) r/.s 



i 



C — IOO 

c+ioo 



J e s s^ ^2 log 2 - 2 log(Vs + x) + log ds 



2m 

C — lOO 

c+ioo 



= 2 log 2 + ^(7) - / e s s^\og{^ + x) ds, (62) 

27T2 J 

c—ioo 

where the last equality follows with aid of the identities 

c+ioo c+ioo 

TO]- 1 = f e ' s ~ 7 ^> and = ^ / e fl S ^log( S )d S , »( 7 ) > 0,c> 0, 

c—ioo c—ioo 

By comparing (|62|) with (|50|) , the formula ( |6l]) follows. There is independent confirmation 
for this result and consequently an alternative proof of the important summation flSUp. Let 
us denote the integral representation in (|62|) by 

c+ioo 

fJx) = / e s s" 7 log(v / s + x) ds (63) 

2m J 

c—ioo 

For this particular representation it is clear that / 7 (0) = If^y- By differentiating ( |63"D with 
respect to x, we get 

c+ioo c+ioo 

i_/ ( x ) = J_ / eV-T-lh--)- 1 ^- A / e's-T-^l--)" 1 ^ 
rfx J ^ v J 2m K s J 2m K s J 



C — lOO C— lOO 

1 „ 1 2 ^ 



1 F 1 (l;- + r ,x 2 )-—-— 1 F 1 (l;l + r ,x z ), 7 > -1 (64) 



r( 7 +|r 1V '2 " ' r( 7 + i) 

where we have used the identity JT 



c+ioo 

ff-i 1 /" A; 

iF 1 (a;w,kt) = — / e s *s^(l - -)~ a ds, u>0,s>k 
2m J s 



C— lOO 
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By integrating (|6^) , using 

J 1^(1; i + 7; t 2 ) = a; 2 F 2 (1, ^; ^ + 7; x 2 ) 



and 



/*x 2 
y txFxil; 1 + 7; t 2 ) dt = X — 2 F 2 (1, 1; 2, 1 + 7; x 2 



we get 



f J x ) = + x 1x 2 F 2 (l, -; -, - + 7; x 2 ) 7 ^-— 2 F 2 (l, 1; 2, 1 + 7; x 2 ) 

1/71 ; 2r( 7 ) 21X7 + !) 2' 2' 2 ; 2r( 7 + i) v ' ' ' h ' 

which proves (^l|) and consequently confirms (|50D - Finally, we mention two identities, among 
many, that follow easily from our work. It is known || 



L^ 2 F x (-n, a; (3; z)t n = (1 - t) x 2 F 1 (-\, a; (3; --^-] 
n\ 1 — t 



n=0 

we may now add, result of Lemma 5, that 

" (" A )« * f „ „.. /3 .^_ r (A + «)r(^)^ 

rc=0 



E ^ l( _ n , a; ^ H L|™. (65) 



Furthermore, by setting 7 = b in (49) and for < y < 1, we have the following identity 



^44^XM'- (66) 
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7 Conclusion 



This paper reports results in an on-going study of the spectrum of an oscillator perturbed 
by a term of the form X/x a . This problem exhibits particularly interesting features when it 
is "singular", that is to say when a > |. In this regime, standard perturbation theory breaks 
down in the sense that the perturbation is not turned off smoothly as A — > 0. We use as 
an un-perturbed basis problem the soluble case a — 2. Although this is is not singular it 
nevertheless provides a basis that is 'closer' to the original singular problem than would be 
the eigenfunctions, say, of an unperturbed oscillator. The perturbation series generated in 
this way, both for the eigenvalues and for the eigenfunctions, lead us in turn to summation 
problems of special series of hypergeo metric functions. We are happy to report that, by the 
present work, all the difficulties of summing these series have now been overcome for the 
case where a is an integer. Consequently we are able to provide explicit formulas for the 
coefficients in the original perturbation series. 
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